Galaxies 2013, xx, 1-x; doi: 10.3390/- 



OPEN ACCESS 



S 



galaxies 



ISSN 2075-4434 

www.mdpi.com/joumal/galaxies 
Article 

Exact expressions for the pericenter precession caused by some 
Dark Matter distributions and constraints on them from orbital 
motions in the Solar System, in the double pulsar and in the 
Vn Galactic center 

pq Lorenzo lorio ^'* 

>^: 

J2 ■ ^ Italian Ministry of Education, University and Research (M.I.U.R.), Fellow of the Royal Astronomical 



Society (F.R.A.S.), Viale Unita di Italia 68, 70125, Bari (BA), Italy. Tel. +39 3292399167 
Received: xx / Accepted: xx / Published: xx 



^ 



(^- Abstract: We analytically calculate the secular precession of the pericenter of a test particle 

orbiting a central body surrounded by a continuous distribution of Dark Matter (DM) by 
^ ' using some commonly adopted spherically symmetric density profiles for it. We obtain exact 

^ I expressions without resorting to a-priori simplifying assumptions on the orbital geometry 

cn ! of the test particle. Our formulas allow us to put constraints on the parameters of the DM 

• ! distributions considered in several local astronomical and astrophysical scenarios, such as the 

O ■ Sun's planetary system, the double pulsar, and the stellar system around the supermassive 

black hole in Sgr A*, all characterized by a wide variety of orbital configurations. As far 
as our Solar System is concerned, latest determinations of the supplementary perihelion 
K/ • precessions Atu with the EPM2011 ephemerides and the common power-law DM density 

profile pBuir) = po^"^A^ yield 5 x 10^ GeV cm-^ (7 = 0) < po < 8 x 10^ GeV cm-^ (^ = 
4), corresponding to 8.9 x 10~^^ g cm^'' < Po < 1-4 x 10^^*^ g cm""^, at the Saturn's distance. 
From the periastron of the pulsar PSR J0737-3039A and the same power-low DM density, 
one has 1.7 x 10^^ GeV cm~^ (7 = 0) < po < 2 x 10^^ (7 = 4) GeV cm~^, corresponding 
to 3.0 X 10^* g cm~'^ < Po < 3.6 x 10~^ g cm~'^. The perinigricon of the SO-2 star in Sgr 
A* and the power-law DM model give 1.2 x 10^^ GeV cm~^ (7 = 0) < po < 1 x 10^^ (7 = 
4, A = rjnin) GeV cm~^, corresponding to 2.1 x 10~^^ g cm"^ < Po < 1-8 x 10"*^ g cm~^. 
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1. Introduction 



Latest results [1] from the ESA's Planck mission [2] have yielded a significative revision of the 
currently accepted content of non-baryonic Dark Matter (DM) of the Universe, which would now 
amount to about 26.5 ± 1.1% with respect to the previously accepted value of 23.2 ± 1.7% from 
WMAP data. Indeed, the Planck-based cold DM density, normalized to the critical density, is [1] 
n^h^ = 0.1199 ± 0.0027, with h = 0.673 ± 0.012, while the analysis of nine years of WMAP data 
yielded [3] Qch^ = 0.1138 ± 0.0045, with h = 0.70 ± 0.022. 

Moving to a galactic scale, the Large Area Telescope (LAT) on the Fermi Gamma Ray Space 
Telescope spacecraft recently discovered a gamma-ray excess at the Galactic center [4-6] which may be 
due to DM annihilation phenomena. The center of the Milky Way is one of the Galactic regions where 
a high DM density is expected. As a consequence, if DM was made of self-annihilating constituents 
into the Standard Model (SM) particles, signals of DM annihilations should come primarily from the 
Galactic center. The gamma-ray excess signature recorded by the Fermi-LAT instrument [4-6] is in the 
form of a monochromatic gamma-ray line with an energy of about 130 GeV. Later, a second peak at 1 10 
GeV was discovered [7]; the same double peak spectrum was independently observed also in gamma-ray 
excess from nearby galaxy clusters [8]. The occurrence of a double peak is a generic prediction of DM 
annihilation pattem in gauge theories. 

There is a lingering debate concerning the amount of DM in our Solar System [9-12]. Recent, 
controversial estimates by Moni Bidin et al. [10] of the DM local density in the neighbourhood of 
our Solar System, based on certain assumptions about the vertical velocity dispersion of old tracer stars 
of the thick disk at the Sun's neighborhood, point towards a value of 

Pdm < 0.04 GeV cm-l (1) 

The bound of eq. (1) is smaller by about one order of magnitude than the usually accepted estimate 

Pdm ~ 0.3 GeV cm'^ (2) 

The figure in eq. (2) has been considered as still valid in other recent studies [11,12] criticizing Moni 
Bidin et al. [10]. Indeed, Bovy and Tremaine [11], who retained the assumptions by Moni Bidin et al. 
[10] implausible, obtained 

Pdm = 0.3 ± 0. 1 GeV cm"^ (3) 

by using a different hypothesis concerning the mean azimuthal velocity of the stellar tracers, while Nesti 
and Salucci [12] inferred 

Pdm = 0.43 ±0.11 ±0.10 GeV cm"^ (4) 

with a different method, independently of the Galactic DM density profile. An accurate knowledge 
of the DM distribution in the neighbourhood of the Solar System is relevant for the attempts aimed 
to directly detect DM particles in laboratory-based experiments such as, e.g., CDMSI [13], CDMSII 
[15], DAMA/Nal [16] and its successor DAMA/LIBRA [17], XENONIO [18] and ZEPLIN III [19]. 
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Interestingly, the silicon (Si) detectors of the CDMS II experiment has recently revealed three WIMP- 
candidate events at a 3a level. [14]. Moving to the Earth's neighbourhood, hints of DM might have been 
detected [20] by the Alpha Magnetic Spectrometer (AMS) on the International Space Station, which has 
recently measured an anomalous high-energy positron excess in Earth-bound cosmic rays [21]. 

Thus, it has become even more important to devise independent means to effectively gain information 
on the constitution and the distribution of such a hypothesized, elusive ingredient of the natural world. In 
Section 2, we will analytically calculate the secular perturbations to the motion of a test particle orbiting 
a central body surrounded by a continuous DM distribution. We will adopt some spherically symmetric 
local DM density profiles [22] by computing the corresponding corrections t/oM to the Newtonian two- 
body potential through the Poisson equation. Then, we will evaluate them onto the unperturbed Keplerian 
ellipse of the test particle, and we will average them over one orbital period Pb of the test particle in order 
to use the standard Lagrange planetary equations [23] for the variation of the Keplerian orbital elements. 
We will, thus, be able to straightforwardly infer exact expressions of the orbital rates of change due to 
the different DM distributions considered without resorting to any a-priori simplifying assumptions on 
the orbital configuration of the test particle. Such formulas will be useful to constrain the DM local 
distributions in different astrophysical and astronomical scenarios such as planetary systems, binaries 
hosting compact objects, stellar systems around supermassive black holes, etc. characterized by a wide 
variety of orbital geometries. Only more or less approximate expressions for the DM-induced pericenter 
precession exist in literature. 

In Section 3 we will apply them to the planets of the Solar System, to the double pulsar, and to the 
stars at the Galactic center in Sgr A*. 

Section 4 summarizes our findings. 

2. Orbital precessions for various DM density profiles 

So far, several authors have put dynamical constraints on the DM distribution within our Solar 
System from orbital motions of its major bodies [24-36]; for some effects of Solar System's DM on 
the propagation of electromagnetic waves, see [35,37]. 

In many cases, more or less approximate expressions for the anomalous perihelion precession induced 
by certain spherically symmetric DM distributions were used, in particular by considering nearly circular 
orbits. In this Section, we will overcome such a restriction by calculating exact expressions which can, 
thus, yield more accurate constraints in view of the increasing level of accuracy in determining the orbits 
of some of the major bodies of the Solar System. Moreover, our results can be used also with systems 
characterized by highly eccentric orbits such as extrasolar planets, binaries hosting compact objects, and 
stellar systems orbiting supermassive black holes lurking in galactic nuclei. 

2.1. Exponential DM density profile 

By choosing an exponentially decreasing DM profile [36] 

PDM(r) = Poexp (^--j , (5) 
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where po is the density at a distance r equal to the characteristic scale length A, the Poisson equation 
yields 



C7„„(.) = -±l£^A' 



exp 



The density profile of eq. (5) is a particular case of the Einasto profile [38], often adopted to describe 
DM halos in galaxies [22]. 

The effect of eq. (6) on the orbital motion of a test particle moving around a central body of mass 
M surrounded by a DM continuous distribution characterized by the density profile of eq. (5) can be 
computed perturbatively by assuming eq. (6) as a small correction to the Newtonian monopole Uj<i{r) = 
—GM/r. The average of eq. (6) over one orbital revolution of the test particle turns out to be 



(7) 



{Ubm) = _ ^^^^o^ ^^P( a) [3 exp (^) A - (a + 2 A) Jq (x) + aeh {x 
a I \A/ 

In eq. (7), Iq (x) , h (x) are the modified Bessel functions of the first kind, and x = ae/\, where a, e are 
the semimajor axis and the eccentricity of the orbit of the test particle, respectively. 

The Lagrange equation for the variation of the longitude of the pericenter zu [23] and eq. (7) yield the 
following secular precession 

dzu\ ATiGpoXVT^^exp (-f ) [-h (x) + eh (x)] 

) (o) 



dt / aeny, 



where nb = \/GM/a^ = 2ti/ P\, is the Keplerian mean motion of the test particle. Although Pitjev 
and Pitjeva [36] considered the density profile of eq. (5), they did not explicitly calculate the resulting 
pericenter precession. The precession of eq. (8) is an exact result in the sense that neither a-priori 
simplifying assumptions on e nor on A were assumed. Note that eq. (8) is not singular in e since, in the 
limit e — )■ 0, it reduces to 

d^\ ^ 27rGpoexp(-f) ^ ^^^ 

dt I rib 

In the case of an infinite length scale, i.e. for an uniform mass density, eq. (8) reduces to the known 
result [27-29,36] 

Aw\ 2nGpoy/l - e^ 

dt I rib 

2.2. Power-law DM density profile 

By assuming a power-law DM distribution profile 

PDM(r) = Po(^) ,7>0 (11) 

usually used for the galactic halos and in several DM-related studies [22,39], the Poisson equation yields 

TT ( \ 47rGpo 2-7 .nx 

^""^^^= (3-7)(2-7)aV '■ (12) 
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The orbital effects of eq. (11) have been computed more or less explicitly and at various levels of 
approximations in [30,32,34,35,40,41]. Actually, also the averaged pericenter precession induced by 
eq. (12) can be exactly computed with the Lagrange pertubative scheme [23] without resorting to any 
simplifying assumptions on the eccentricity e of the test particle. 

It turns out to be 

du7\ TrGpo v^ / ,x ..ox 

-Jt) = r VPj e;7,A , (13) 

where the cumbersome expressions of the coefficients Pj(e;7, A), j = 1,2,... 16 are explicitly 
displayed in Appendix A. It is worthwhile noticing that eq. (13) is not any sort of somewhat truncated 
series expansion, and that the coefficients in Appendix A are exact in the sense that no a priori simplifying 
assumptions were assumed in calculating them. 

By performing a series expansion of eq. (13) in powers of e one obtains 



/dw\^ nGpo pVfo 1 
\ dt / n^ \a) \ 


'7(7- 1)- 
4 


-4" 


e' 


The first term in eq. (14) yields a perihelion shift per orbit 




/dw\ /27i\ 
Aw = {^) — 
\dt / \nkj 




in agreement with that by Frere et al. [32]. 









\+0 (e^) . (14) 



(15) 



3. Confrontation with the observations 

In this Section, we will use our analytical predictions of eq. (8) and eq. (13) to infer constraints on po 
as functions of A and 7 in some astronomical and astrophysical scenarios. 

3.1. Planets of the Solar System 

By following the approach adopted in several researches on DM [26-36] and on other putative non- 
standard dynamical effects [27,42] in our Solar System, we will confront our theoretically predicted 
pericenter precessions of eq. (8)andeq. (13) with the latest determinations of the admissible ranges Azu 
for any possible anomalous perihelion precessions obtained by fitting up-to-date dynamical models to a 
centennial record of planetary observations of several types. The term "anomalous" refers to the standard 
Newtonian-Einsteinian dynamics, fully modeled in the most recent planetary dynamical theories. At 
present, two independent teams of astronomers are engaged in producing, among other things, such 
corrections Aw to the planetary perihelion precessions [36,43]; here we will adopt the latest estimates 
by Pitjev and Pitjeva [36] , summarized in Table 1 , which are based on the most recent version of the EM? 
ephemerides. DM was not explicitly modeled in the EPM routines; thus, the supplementary precessions 
of Table 1 are well suited to put constraints on it in a phenomenological, model-independent way. Here 
and in the rest of the paper, "model-independent" has to be intended as independent of the nature and 
of the physical properties of the DM particles. On the other hand, the constraints inferred in this paper 
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Table 1. Supplementary precessions Aw of the longitudes of the perihelion for some planets 
of the Solar System estimated by Pitjev and Pitjeva [36] with the EPM2011 ephemerides. 
They processed 676804 data of various kinds covering about one century (1913-2010) by 
estimating more than 260 parameters. While some tracking data from Cassini were used, 
those from Messenger, currently orbiting Mercury, were not included in the data record. 
Pitjev and Pitjeva [36] fully modeled all standard Newtonian-Einsteinian dynamics, apart 
from the general relativistic Lense-Thirring effect caused by the Sun's rotation. However, its 
expected magnitude from general relativity is smaller than the uncertainties quoted here for 
the planets used in the text (Earth, Mars, Saturn), so that it does not impact our results. . The 
units are milliarcseconds per century (mas cty~^). The errors released a^xb^ in nias cty^^ as 
well, generally exceed the formal, statistical ones by several times [36]. 



Planet 


Atxj (mas cty ^) 


cTAtii (mas cty ^) 


Mercury 


-2.0 


3.0 


Venus 


2.6 


1.6 


Earth 


0.19 


0.19 


Mars 


-0.020 


0.037 


Jupiter 


58.7 


28.3 


Saturn 


-0.32 


0.47 



depend on how DM clusters in space, i.e. on the specific density profiles adopted. Moreover, by assuming 
that the whole ranges Aw in Table 1 are entirely due to DM, the resulting constraints will be relatively 
"generous", mitigating the risk of getting unrealistically tight bounds. A complementary approach which 
could be followed, recommended by some researchers [44] when non-standard effects such as DM are 
involved, consists in explicitly modeling DM and fitting such ad-hoc modified dynamical theories to the 
same data set in order to estimate a dedicated solved-for parameter in a least-square way. One could 
try with different DM density profiles to establish the one providing the best fit to the data. Such a 
strategy would be particularly meaningful in the case where non-zero extra-precessions, at a statistically 
significant level, should be extracted from future observations. 

3.1.1. Constraints for the exponential density profile 

Here, we deal with the exponential density profile of eq. (5). 

In Figure 1 we plot the upper bounds on po, in GeV cm"'^, as a function of the characteristic scale 
length A over the extension of the orbits of the Earth, Mars and Saturn. We compare eq. (8) to the figures 
in Table 1. The tightest bounds come from the perihelion of Saturn; 1.7 x 10^ GeV cm~'^ < Po ^ 
1.9 X lO'' GeV cm"^ corresponding to 3.0 x 10"^° g cm'^ < Po < 3.4 x 10"^° g cm'^. 

3.1.2. Constraints for the power-law density profile 

Here, the power-law distribution of eq. (11) is considered. 
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Figure 1. Upper bounds, in GeV cm^'^ (1 Gev cm^'^ is equivalent to 1.78 x 10"^^ g cm^'^), on 
the DM density parameter po as a function of the characteristic scale length A of eq. (5). The 
supplementary perihelion precessions of Table 1 for Earth, Mars and Saturn were adopted 
along with the analytical prediction of eq. (8). The characteristic length A, in AU, is assumed 
to vary from fmin = a(l — e) to r^ax = a(l + e). 
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Figure 2 displays the bounds on po, in GeV cm~'^, as a function of the characteristic scale length A 
and of the parameter 7 for the Earth, Mars and Saturn. Also in this case, A covers the orbit extensions 
of the planets considered, while [32,35] < 7 < 4. Also in this case, Saturn yields the tightest bounds: 
5 X 10^ GeV cm~^ (7 = 0) < po < 8 x 10^ GeV cm~^ (7 = 4) corresponding to 8.9 x 10~^^ g cm~^ < 
Po < 1-4 X 10~^° g cm~'^. They are about one order of magnitude better than those inferred in Section 
3.1.1. . 

3.2. The double pulsar 

Thanks to their high matter density, astrophysical compact objects such as neutron stars should be able 
to efficiently capture significative amounts of DM [45-55] depending on the nature of its particles, on 
how it clusters in space and on the history of the stars themselves. Concentrations of DM around neutron 
stars may give rise to phenomena such as gamma-ray emission [54,55] which, on the one hand, depend 
on the physical properties of DM particles themselves [56] and, on the other hand, can be generated also 
by other physical mechanisms [57,58]. Thus, it is valuable to constrain the DM density around binaries 
hosting neutron stars in a phenomenological and model-independent way. 

To this aim, we look at the double pulsar PSR J0737-3039A/B [59,60] . As in Section 3. 1 , we consider 
both eq. (5) and eq. (11). In this case, we use the periastron precession of PSR J0737-3039A, which 
is nowadays determined with an accuracy of 6.8 x 10~^ deg yr"^ from timing measurements [61]. 
Nonetheless, care is required in straightforwardly using such a figure in comparisons with theoretical 
predictions of non-standard effects to infer constraints on them. Indeed, a larger uncertainty should, 
actually, be considered on the periastron rate of PSR J0737-3039A because of the mismodeling in its 
IPN periastron precession [62]; it may be as large as 0.03 deg yr~^ [62]. In Figure 3 the bounds for both 
eq. (5) and eq. (1 1) are shown. The power-law profile provides the tightest constraints. Indeed, over the 
extension of the double pulsar system, it yields 1.7 x 10^^ GeV cm"^ (7 = 0) < po < 2 x 10^^ (7 = 4) 
GeV cm~'^ corresponding to 3.0 x 10^^ g cm"'^ < Po < 3.6 x 10^^ g cm~'^, while the bounds from the 
exponential density profile are of the order of 5 x 10^^ GeV cm"'^ corresponding to 9 x 10^^ g cm~^. 

3.3. The stellar system orbiting the Galactic black hole in SgrA* 

The issue of the DM distribution at the center of galaxies in presence of supermassive black holes 
has been treated in, e.g., [40,41,63-66]. It depends on several factors such as [63] the galactic halo 
density profile itself [67] and the nature of the DM particles themselves. Moreover, as pointed out in 
[66], significant differences in the final DM distribution close to the black hole are found depending on 
the theoretical scheme adopted for the calculation. 

So far, the presence of DM at the center of the Milky Way has been indirectly guessed from 
certain phenomena such as gamma-ray emissions interpreted as DM annihilation [68,69]. However, 
such interpretations are not free from more or less tested assumptions about several other concurring 
physical phenomena [70,71]. Recently, the Fermi satellite detected a gamma-ray excess at the Galactic 
center which may be due to DM annihilation phenomena [4-6]. Thus, it is important to constrain DM 
distributions at the Galactic center in a dynamical, model- independent way. 
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Figure 2. Upper bounds, in GeV cm^'^ (1 Gev cm~^ is equivalent to 1.78 x 10"^"^ g cm~^), 
on the DM density parameter po ^s a function of the parameters A and 7 of eq. (11). The 
supplementary perihelion precessions of Table 1 for Earth, Mars and Saturn were adopted 
along with the analytical prediction of eq. (13). The characteristic length A, in AU, is 
assumed to vary from Tmin = a(l — e) to r^ax = (^(l + e). 
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Figure 3. Upper bounds, in GeV cm^'^ (1 Gev cm^'^ is equivalent to 1.78 x 10^ g cm^'^), 
on the DM density parameter po in the PSR J0737-3039A/B system as a function of A (top 
panel) of eq. (5), and of A and 7 of eq. (11) (bottom panel). The characteristic length A, in 
AU, is assumed to vary from r^i^ = a{l — e) to rmax = o-i^ + e). A conservative periastron 
rate uncertainty of 0.03 deg yr~^ [62] was assumed for PSR J0737-3039A. It was compared 
to the theoretical predictions of eq. (8) (top panel) and of eq. (13) (bottom panel). 
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Figure 4. Upper bounds, in GeV cm^'^ (1 Gev cm~'^ is equivalent to 1.78 x 10~ g cm~'^), 
on the DM density parameter po in the Sgr A*— SO-2 system as a function of A (top panel) 
of eq. (5), and of A and 7 of eq. (11) (bottom panel). The characteristic length A, in AU, 
is assumed to vary from rmin = a(l — e) to rmax = a(l + e). A conservative perinigricon 
rate uncertainty of 0.6 deg yr~^ was assumed for SO-2. It was compared to the theoretical 
predictions of eq. (8) (top panel) and of eq. (13) (bottom panel). 
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To this aim, we will look at the orbital motions of the stellar system [39,72,73] revolving about the 
supermassive black hole located at the Galactic center in Sgr A* [74]. So far, available data cover 
one full orbital revolution for two stars: SO-2, having an orbital period P^ = 16.17 ± 0.15 yr and 
eccentricity e = 0.898 ± 0.005 [73], and SO-102, characterized by Pb = 11-5 ± 0.3 yr and eccentricity 
e = 0.68 ± 0.02 [73]. The mass of the black hole is M, = (4.1 ± 0.4) x IO^Mq [73]. From such 
figures, a naive, order-of-magnitude evaluation on the accuracy which could be reached in determining 
the stellar perinigricon precessions can be made; for SO-2 we have a^j ~ 0.6 deg yr~\ while for SO-102 
it can be inferred o-^j ~ 4.8 deg yr~^. Then, we will adopt SO-2 as probe to constrain DM density as 
done in the previous Sections. In view of its high eccentricity, the exact formulas of eq. (8) and eq. 
(13) are well suited for such a task. From Figure 4 it turns out that the bounds for the exponential 
density vary sensibly over the SO-2 orbit extension in view of its high eccentricity; indeed, it is po ~ 



4 X 10 GeV cm •^, corresponding to 7 x 10 g cm '^, for A 



while it reduces to po 



2 X 
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10^^ GeV cm^'^, corresponding to 3 x 10^^^ g cm"^, for A = rmax- As pointed out in [39], as far as 
supposedly baryonic DM is concerned, both theoretical [75,76] and observational [77] motivations for 
considering a power-law profile such as eq. (11) at the galactic centers exist in literature. The same 
kind of potential was considered for non-baryonic DM as well [41,64]. According to Figure 4, it yields 
1.2 X 10^3 GeVcm-3 (7 = 0) < po < 1 x 10^^ (7 = 4, A = rmin) GeV cm'^ corresponding to 
2.1 X 10"^^ g cm~^ < Po < 1-8 X 10"^ g cm~^ 

For an earlier idea of constraining the DM concentration near the Galactic center from the perinigricon 
ofSO-2, see [40,41]. 

4. Summary and conclusions 

Recent results concerning the amount of non-baryonic DM, both at global and local scales, stress the 
need of further deepening the research of accurate and independent strategies to gain information about 
the distribution of such a hypothetical key ingredient of the natural world. 

To this aim, we looked at the effects induced by some spherically symmetric DM density profiles on 
the motion of a test particle orbiting a localized body surrounded by an extended DM distribution. In 
view of the increasing accuracy in the determination of the orbits of the planets of our Solar System and 
of the possibility of looking also at different astronomical and astrophysical laboratories characterized by 
a wide variety of orbital configurations, we analytically calculated the DM-induced pericenter precession 
without resorting to any a-priori simplifying assumptions concerning the orbital geometry of the test 
particle. In this respect, our results, obtained perturbatively by means of the Lagrange planetary 
equations, are exact, being valid for any value of the eccentricity of the orbit of the test particle. 

We considered an exponentially decreasing profile pBuir) = po exp (— r/A) and a standard power- 
law model pBuir) = p^r^'^X' . We compared our analytical predictions with the latest observational 
determinations for some planets of our Solar System obtained with the EPM2011 ephemerides, the 
double pulsar, and the SO-2 star orbiting the supermassive black hole in Sgr A*. The tightest constraints, 
obtained for the power-law model, came from the supplementary precessions of the planetary perihelia. 
We obtained 5 x 10^ GeV cm~^ (7 = 0) < po < 8 x 10^ GeV cm~^ (7 = 4), corresponding to 
8.9 X 10"^^ g cm"'^ < Po < 1-4 x 10~^° g cm"'^, at the Saturn's distance. From the periastron of 
PSR J0737-3039A we inferred 1.7 x 10^^ GeV cm'^ (7 = 0) < po < 2 x 10^^ (7 = 4) GeV cm-3, 
corresponding to 3.0 x 10"*^ g cm~^ < Po ^ 3.6 x 10~^ g cm~^. The perinigricon of the SO-2 star in Sgr 
A* gave 1.2 x 10^^ GeV cm"^ (7 = 0) < po < 1 x 10^'' (7 = 4, A = rmin) GeV cm"^, corresponding 
to 2.1 X 10~i^ g cm-3 < Po < 1.8 x 10"^ g cm-^ 

Our results can be used in future when new, more accurate data will be collected and processed. As a 
complementary approach which could be followed, DM dynamical effects should be explicitly modeled 
in the softwares used to reduce the planetary observations, and dedicated solved-for parameters should 
be estimated in fitting the newly constructed models of the forces acting on the planets to given data 
records. 

Appendix 

A. The coefficients of the pericenter precession due to the power-law density Below, the exact 
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expressions of the coefficients p^ (e, ; 7, A) , j = 1, 2, . . . , 16 of the pericenter precession of eq. (13), 
induced by eq. (1 l)-eq. (12), are listed. No simplifying assumptions on the eccentricity e were assumed 
in computing them. They contain the hypergeometric function 2F1 (g, w; b; d) [78]. 



2(l + e)-^A^ 
Pi ^ 



P2 



15 (-1 + ey eyr^e2 (6 - 57 + 72)^ 
(7 - 3)(7(7(167(87 - 61) + 2687) - 3204) + 1405)e^- 

- (7 - 2)(7(167(7(87 - 57) + 140) - 2265) + 815)e^+ 
+ (7(7(7(4647 - 3917) + 11838) - 15306) + 7155)e^+ 
+ 2(7(7(247(57 - 38) + 2813) - 4095) + 2230)6^+ 

+ (7((5099 - 9157)7 - 8991) + 5095)6^+ 
+ 15(7(357 -151) + 154)e- 

- 5(37((7 - 8)7 + 24) - 73) 



15 (-1 + e) eVr^e2 (6 - 57 + 72)^ 
2(7 - 3)(7 - 2)(7(27(87(87 - 55) + 1093) - 2421) + 985)e^- 

- (7 - 2) (7(47(7(87(87 - 73) + 2065) - 3593) + 12571) - 4325)e*^+ 
+ (7(7(27(127(527 - 485) + 21371) - 77699) + 69856) - 24575)e^+ 
+ (7(87(57(37(47 - 47) + 616) - 6337) + 49979) - 19385)e^- 

- 2(7(7(7(9907 - 7771) + 21801) - 26063) + 11335)6^+ 
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+ (7(7(60(35 - 37)7 - 8891) + 15549) - 9320)e^+ 
+ 5(37 - 5)(7(307 - 119) + 95)e- 
- 5(97(27 - 9) + 89) 

}.F,(1,7;2;^), (A.2) 



^ 27(l + e)"^A^ 

^^~ 5(-l + e)v^r^^(6-57 + 72)2^ 



(47(7(87 - 53) + 104) - 225)e=^ 
+ 2(47-ll)(47-5)e2- 
- 5(47(37 - 14) + 63)e+ 
+ 10 



} 2Fi(^,1+7;2;t^), (A.3) 



Pi 



2' " '1 + e^ 

7(l + e)^n^ 



5 (-1 + ey Vr^e2 (6 - 57 + 72)^ 
(7(7(87(87 - 67) + 1669) - 2256) + 1070)e^+ 
+ (7(47(47 - 45) + 489) - 335)e^- 

- (27 - 3) (37 - 10)(207 - 47)6^- 

- 5(7(127- 55) + 65)e+ 
+ 5(37-8)(37-7) 
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} 2Fi(^,l+7;3;^), (A.4) 



P5 = 



5 (-1 + ey eVT^^ (6 - 57 + 7^)^ 
(7 - 3) (47(7(87 - 53) + 104) - 225)e^+ 
+ (7(1839 - 47(7(87 - 85) + 311)) - 895)e^+ 
+ (7(47(97 - 71) + 575) - 255)6^+ 
+ (27 - 5)(7(307 - 107) + 123)6^- 
- 10(7(127- 55) + 60)e+ 
+ 10 



}2i^i(-^,7;l;^), (A.5) 



2(l + e)'-n 

P6 = 



5 (-1 + e)^ evT^^e^ (6 - 57 + 7^)^ 
(7 - 3)(7(7(87(87 - 67) + 1669) - 2256) + 1070)e^+ 
+ (7(7(7(8(101 - 87)7 - 4025) + 9781) - 11407) + 4950)e^+ 
+ (7(7(27(927 - 727) + 4311) - 5616) + 2600)e^+ 
+ (7 (27 (607^ - 6697 + 2483) - 763l) + 4300) e^+ 
+ (7((1891 - 2557)7 - 4409) + 3210)e^- 
- 5(37(7(37 - 19) + 42) - 94)e+ 
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+ 5(37-8)(37-7) 



}2i^i(-^,7;2;^), (A.6) 



27(l + e)-"A^ 
^^ 15(l-e2)^/'e(6-57 + 72)2^ 



(7(7(167(87 - 61) + 2687) - 3204) + 1405)e^ 
+ (7((447 - 647)7 - 939) + 590)e^- 
- 2(7(7(1207 - 649) + 1066) - 525)6^+ 
+ 10(37 -8)(77-ll)e2+ 
+ 5(3(7- 8)7 + 37)e+ 
+ 15(7-3)(7-2) 



}2Fi('|l+7;2;^), (A.7) 



7(l + e)"'^''^A^ 
^' 15(l-e)'/'(6-57 + 72)' 



256e^7^- 

- 32e^(e(71e + 6) + 15)7^+ 

+ 4e(e(e(e(1973e + 348) + 870) + 60) + 45)7^- 

- 2(e(e(e(e(6793e + 1793) + 4446) + 826) + 465) + 45)7^+ 
+ (e(e(e(e(11654e + 3925) + 9608) + 3518) + 1370) + 405)7- 
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5(e(e(e(e(788e + 317) + 764) + 458) + 104) + 89) 



} 2Fi(|l+7;3;^), (A.8) 



2(l-e)-^A^ 

P9 - 



15 (1 + ef^^ eVT^^ (6 - 57 + -f^f 
(7 - 3)(7(7(167(87 - 61) + 2687) - 3204) + 1405)e^+ 
+ (7 - 2)(7(167(7(87 - 57) + 140) - 2265) + 815)e^+ 
+ (7(7(7(4647 - 3917) + 11838) - 15306) + 7155)e^- 

- 2(7(7(247(57 - 38) + 2813) - 4095) + 2230)6=^+ 
+ (7((5099 - 9157)7 - 8991) + 5095)6^- 

- 15(7(357 -151) + 154)e- 

- 5(37((7 - 8)7 + 24) - 73) 



}2i^i(^,7;i;3^), (A.9) 



2(l-e)-^A^ 



15 (1 + e)^/^ e^/T^^ (6 - 57 + 72)^ 
2(7 - 3)(7 - 2)(7(27(87(87 - 55) + 1093) - 2421) + 985)e^+ 
+ (7 - 2) (7(47(7(87(87 - 73) + 2065) - 3593) + 12571) - 4325)6*^+ 
+ (7(7(27(127(527 - 485) + 21371) - 77699) + 69856) - 24575)e^+ 
+ (7(-87(57(37(47 - 47) + 616) - 6337) - 49979) + 19385)e^- 
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- 2(7(7(7(9907 - 7771) + 21801) - 26063) + 11335)e^+ 
+ (7(7(607(37 - 35) + 8891) - 15549) + 9320)6^+ 
+ 5(37 - 5)(7(307 - 119) + 95)e+ 
+5(97(27 - 9) + 89) 

}.F,(1,7;2;-^), (A.10) 

27(l-e)-'/'-^A^ 



P12 



5(1 + 6)^/^(6-57 + 72)^ 
32e^7^- 

- 4e(e(53e + 8) + 15)7^+ 
+ 8e(4e(13e + 4) + 35)7- 

- 5(e(e(45e + 22) + 63) + 2) 

}.F,(l,l+7;2;3^), (A.11) 

7(l-e)-n^ r 



5(1 + ef^^ ^JV^e (6 - 57 + 72)^ 
(7(7(87(87 - 67) + 1669) - 2256) + 1070)e^+ 
+ (7(4(45 - 47)7 - 489) + 335)e^- 
- (27 - 3) (37 - 10)(207 - 47)6^+ 
+ 5(7(127- 55) + 65)e+ 
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+ 5(37-8)(37-7) 



} 2Fi(i,l+7;3;-^j, (A.12) 



^'' ~ ~5e(l + e)'/'(6-57 + 72)' 



(7 - 3) (47(7(87 - 53) + 104) - 225)e^+ 
+ (7(47(7(87 - 85) + 311) - 1839) + 895)e^+ 
+ (7(47(97 - 71) + 575) - 255)e=^+ 
+ (7(4(91 - 157)7 - 781) + 615)6^- 

- 10(7(127- 55) + 60)e 

- 10 



}2i^i(-^,7;i;3^), (A.13) 



2(-l + e)'(l-e)"^A^ 
^'^~ 5e(l-e2)^/2(6-57 + 72)2^ 



(7 - 3) (7(7(87(87 - 67) + 1669) - 2256) + 1070)e^+ 

+ (7(7(7(87(87 - 101) + 4025) - 9781) + 11407) - 4950)e^+ 

+ (7(7(27(927 - 727) + 4311) - 5616) + 2600)e^+ 

+ (7 (7631 - 27 (6O72 - 6697 + 2483)) - 4300) e^+ 

+ (7((1891 - 2557)7 - 4409) + 3210)6^+ 
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+ 5(37(7(37 - 19) + 42) - 94)e+ 
+ 5(37-8)(37-7) 



27(l-e)~n^ 
15e(l-e2)^/2(6-57 + 72)' 






(7(7(167(87 - 61) + 2687) - 3204) + 1405)e^ 
+ (7(7(647 - 447) + 939) - 590)e^- 

- 2(7(7(1207 - 649) + 1066) - 525)6=^- 

- 10(37 -8)(77-ll)e2+ 
+ 5(3(7-8)7 + 37)6- 
-15(7-3)(7-2) 



'3 2e 

} 2Fi(-,l+7;2;3Y^), (A.15) 



7(l-e)~^A^ 
15 (1 - e)'/' (1 + ef^ (6 - 57 + 72)^ 



^^^ ^3/2 n ^ .^5/2 /R _ c-^, ^ ^ 0^2 ^ 



2(7 - 2)(7(27(87(87 - 55) + 1093) - 2421) + 985)e^+ 
+ (7(27(247(47 - 29) + 1793) - 3925) + 1585)e^- 
- 4(7(37(107(47 - 29) + 741) - 2402) + 955)e^+ 
+ (2290 - 27(27(607 - 413) + 1759))e2+ 
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+ 10(7(37(67 - 31) + 137) - 52)e+ 
+ 5(97(27 - 9) + 89) 

}.f.(|.i + o;3i-i^). (A.16) 
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